Asymmetric localization in disordered Landau bands 
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We show that due to the Landau band mixing the eigenstate localization within the disordered 
bands get an asymmetric structure: the degree of localization increases in the lower part of the band 
and decreases in the upper one. The calculation is performed for a 2D lattice with the Anderson 
disorder potential and we prove that this effect is related to the upper shift of the extended states 
within the band and is enhanced by the disorder strength. The asymmetric localization and the 
energy shift dissapear when the interband coupling is switched off. 
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The localization effect in Landau bands attracted much 
interest since the discovery of the QHE 0, H 0, 0, U B 
IE 13 When the impurity potential is present, the 
initially degenerate Landau levels of the 2D system are 
turning into broad bands. In contrast to the zero mag- 
netic field case when there are no extended states in two 
dimensions 0|, extended states are present at the cen- 
ter of each Landau energy band lj. The early studies 
performed for the continuous Hamiltonian model with 
neglected inter-Landau band mixing show that 

the generic picture of the 2D Landau bands contains lo- 
calized states in the band tails and extended states in the 
middle. However, in a real 2D system the band separa- 
tion may become smaller than the band width so that 
the above approximation is not always valid. It has been 
shown that the mixing between different Landau bands, 
that generally comes from the disorder presence, is re- 
lated to the energy shift of the extended states from the 
central position of the band [H[ll[li[llHH[l3or 
it may have a derealization effect when the states with 
opposite chirality are coupled [Tsl ]. 

The lattice model captures this relevant feature, and 
this happens because the discrete Landau model (initially 
solved by Hofstadter for the pure case 0) takes auto- 
matically into account the inter-band coupling. The en- 
ergy shift between the position of the extended states 
and the peak of the density of states was explicitely cal- 
culated in the lattice model by 22]. Furthermore, for 
the same model, one notes the asymmetric behavior of 
the eigenstates localization in Landau bands. This result 
was reported in 20] by the calculation of the localization 
length within the outer Landau band, and is consistent 
with the asymmetric behavior of the inverse participa- 
tion number within the Landau bands reported recently 
in pH . In this paper we study the origin of this asym- 
metry. By the projection of the 2D discrete Hamiltonian 
on the nondisordered Landau levels we show here that 
this asymmetric behavior of the eigenstates localization 
are due to the interband coupling terms. It will be also 
shown that this effect is intimately related to the shift of 
the extended states. It means that both of the two effects 
disappear when the interband coupling is switched off. 

We use for this the following 2D Landau Hamiltonian 



written in the discrete basis of a 2D rectangular lattice: 
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where V is the Anderson disorder potential: 
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The discrete points (n, m) define the 2D rectangular 
lattice with surface L 2 = N x M and lattice constant a. 
{\n, to}} with n = 1 • • • N and m = 1 • • • M is the discrete 
vector basis and generates the Hilbert space of the one- 
electron states. The periodical boundary conditions are 
used, meaning that \n, M + 1} = \n, 1} and \N + 1, to) = 
|1,to) (2D toroidal geometry). The perpendicular mag- 
netic field in the Landau gauge A = (—By, 0,0) is intro- 
duced by the Peierls substitution in the hopping elements 
along the x — na direction, t — > t exp{27rim0}, where <j) 
is the magnetic flux through the unit cell a 2 of the lat- 
tice measured in the quantum flux units cf>o — h/e. In 
(0) the energies e nm represent the random variables uni- 
formly distributed in the energy interval [— W, W]. W 
is the amplitude of the Anderson disorder potential (or 
disorder strength); t is the energy unit and is set to 1. 
For the commensurate values of the magnetic flux, the 
eigenstate spectrum of the pure system (V=0) exhibits 
the well-known Hofstadter butterfly structure [n| . In the 
numerical calculation we set the flux value as the ratio 
4> = l/p and the system size as L 2 — (integer ■ p) 2 . In 
this case the eigenstates of the 2D system are grouped in 
p discrete Landau bands, every band having a number of 
rib = L 2 /p degenerated eigenstates. To each energy level 
e° , with the band index a = 1 • • • p, correspond rib degen- 
erate eigenvectors j) = \ai) with i = 1 • • • rib- As the 
system has the electron-hole symmetry we concentrate 
on the lower half of the spectrum. 

In the presence of the disorder potential the degenerate 
energy level e° a turns into the broad energy band {e a ;} 
with i = 1 • • • rife. We study the degree of localization of 
the nondegenarate eigenstates \^ a i) by the calculation of 
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the inverse participation numbers IPN which is defined 
as: 

IPN = P ai =J2\ (n, n# ai> | 4 . (3) 

P a i varies from 1/L 2 for the extended states, when the 
electron wave function spreads out over the whole surface 
of the plaquette to 1 for the strong localized states. The 
nature of the eigenstates can also be checked by the cal- 
culation of the variance of the level spacing distribution, 
St [13. 

In this work we put into evidence the role of the inter- 
band interaction. To this end we write the Hamiltonian 
(JTJ in the vector basis of the 2D pure system, {|cu)}: 

P n b p n b 

H (<f>) = EE e °l ai >H + E E ^w|m>(aj| 

a— 1 z—l a — 1 i,j=l 

P n b 

+c J2 H v cri,Pi\oi)(fi3\- ( 4 ) 

In this representation the disorder potential V becomes a 
sum of two terms, corresponding to the inter- and intra- 
band coupling (the second and the third term in re- 
spectively). V a i.fi j a- re the matrix elements of the Ander- 
son potential V written in the basis of the eigenfunctions 
of the ordered system {|m)}. They are random variables 
as well, and their values are proportional to the disorder 
strength W. The coupling constant c is introduced for 
convenience. 

We shall discuss first the situation when the parameter 
c is set to zero in (4) , meaning that only the intraband 
coupling is taken into account. A picture of the disor- 
dered bands for this case is given in Fig. 1. The density 
of states (DOS) for the first two bands is shown in Fig. la 
and its profile has a semielliptic shape. Let Eb be the en- 
ergy where DOS reaches its maximum, which in this case 
is located in the middle of the band. The level spacing 
distribution is calculated by averaging over different dis- 
order configurations in the manner described in p3| . It 
is known that the extended states belong to the unitary 
Wigner-Dyson ensemble (3 = 2 with the variance of the 
level spacing St = 0.42, while the localized ones are dis- 
tributed according to the Poisson law with the variance 
St = 1 24] . The calculation of the level spacing variance 
in Fig. lb shows the presence of extended eigenstates in 
the middle of the disordered bands, for energies E around 
Eb- For energies towards the band edges the states be- 
come localized and there is a continuous crossover from 
unitary Wigner-Dyson distribution to the Poisson distri- 
bution as in [23 . Obviously, the most localized states are 
at the edges of the band, where St increases to higher 
values. We note that in the thermodynamic limit the 
extended states in the band center collapse into a single 
energy level 0- 

We complete the picture of the eigenstates localization 
for uncoupled bands (c = 0) showing in Fig. lc (the red 
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FIG. 1: The density of states (a), the level spacing variance St 
(b), and the inverse participation number (c) for the first two 
Landau bands vs. energy. (I? = 20 2 , W = 2, 0=0.1). DOS 
and St are presented only for the symmetric case with inter- 
band coupling constant c=0. IPN is presented for c= 0, 1 and 

2. In (b), when St = ^ — 1 ~ 0.42, the eigenstates cor- 
respond to the extended states of the unitary Wigner Dyson 
ensemble ((3 = 2). (The configuration average is performed 
over 5000 samples). 

curve) the values of the IPN. By varying the energy from 
the central position of every band, the IPN increases, 
indicating an increased localization of the eigenstates. 
Let E c be the energy with the lowest IPN value in the 
middle of the band. This is the energy of the most ex- 
tended state, where, in the thermodynamic limit, the 
localization-delocalization transition takes place. For the 
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FIG. 2: (a) The inverse participation number P — P m i n and 
(b) the density of states for the first Landau band plotted 
vs. energy at different values of the interband coupling con- 
stant: c=0 (crosses), c=l (stars) and c~2 (dots). (L 2 = 20 2 , 
W = 1, (j> = 0.1). Pmin is the minimum value of the inverse 
participation number, while E c is the energy where this value 
is reached, i.e. P m in = P{E C ). 



case discussed here, when c=0, the critical energy corre- 
sponds to the maximum of the DOS (E c — Ef,), and the 
IPN is symmetric within the band. This is what we call 
the symmetric case. These properties are not preserved 
anymore when c ^ as we could already notice in Fig. lc 
(see the IPN curves for c=l and 2). 

In what follows, we are interested in finding out how 
the localization properties evolve with the interband cou- 
pling c. In Fig. 2 we show the result of the numerical cal- 
culation for the inverse participation number and density 
of states as function of E — E c for three values of the in- 
terband coupling constant c= 0, 1 and 2. For the first 
band, the IPN curves are depicted in Fig. 2a. By the 
definition of E c , the IPN takes the minimum value at 
E — E c = 0. One remarks that the symmetry of the IPN 
is lost for nonvanishing coupling constant c=l and 2, i.e. 
in the case of the band mixing. Compared to the uncou- 
pled case (c=0), the increased values of IPN for E < E c 
at c ^ indicate an increased degree of localization in 
the lower part of the band. The opposite is true in the 
upper part where the states become less localized. 

In Fig. 2b we depict the DOS of the first Landau band. 



For c=l, the band is shifted downwards meaning that 
the maximum of DOS does not correspond to E c , but it 
occurs at a lower energy E{, < E c . The shift increases for 
c=2. It means that the critical energy is moving up in the 
band when the interband coupling is present. This asym- 
metric behavior is preserved for all the bands contained 
in the lower half of the spectrum. 

Once we have established the relation between the 
interband coupling and the asymmetry of the Landau 
bands, we are interested now to learn how this property 
depends on disorder. We keep fixed c=l in Q and in- 
crease the amplitude of disorder W. Since the interband 
coupling in the discrete Hamiltonian is due to the pres- 
ence of the disorder, we expect the shift of the critical 
energy be also dependent on the disorder amplitude. 

Fig. 3 gives the inverse participation number and the 
density of states as function of [E — E c )/W for different 
values of W . One notices in Fig. 3a that, for any disorder, 
in the domain of extended states around E c the inverse 
participation number P(E) can be expressed as P(E) — 
P(E C ) = f((E — E c )/W). Deviations from this law are 
noticed at the band edges. By increasing the disorder 
amplitude W , the lower energy states of the band become 
more localized, but the higher energy states of the band 
become less localized. At the same time, the extended 
states are moving towards the upper edge of the band, 
so that the critical energy E c does not correspond any 
more to the maximum of the density of states. This can 
be seen in Fig. 3b, where one notices the band shift with 
increasing disorder strength. The critical energy shift 
E c — Ef, as function of disorder is shown in Fig. 4 for 
the first two bands in the spectrum. For the present 
model, (E c — Eb)/W exhibits a linear dependence on the 
disorder strength W . We note no significant difference 
of the critical energy shift for the two bands depicted. 
The numerical calculations have been repeated for an 
increased system size, the results being the same. 

In conclusion, we have shown that the band mixing 
gives rise to the asymmetry of the localization properties 
in the Landau bands. The inverse participation number 
(which measures the degree of localization) becomes an 
asymmetric function within the band, indicating that the 
degree of localization increases for the states in the lower 
energy part of the band and decreases for the states in 
the upper part. At the same time the critical energy 
(the most extended state) in each Landau band does not 
correspond to the maximum of the density of states but 
it is shifted to higher energies. These properties are spe- 
cific to the many-band model and are quite different from 
the properties of the one-band model which exhibits only 
symmetrical features. In a large range of energies the 
inverse participation number P{E) — P(E C ) can be ex- 
pressed as a universal function of (E — E c )/W. 

The shift of the extended states from the central posi- 
tion of the Landau band supports the result reported re- 
cently by Shlimak et al [2f| consisting in the displacement 
of the integer values of the filling factor from the middle 
point of the QH plateaus in disordered Si-MOSFET sam- 
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FIG. 3: (a) The inverse participation number P — P m in, 
and (b) the density of states DOS for the first Landau band 
vs. energy at different disorder strength W=0.5 (crosses), 
W=l (stars), W=1.5 (squares), and W=2.0 (dots). (L 2 = 
30 2 , 4> — 0.1, c = 1). Pmin is the minimum value of the inverse 
participation number, while E c is the energy where this value 
is reached, i.e. Pmin = P(E C ). Note that the energies are 
scaled by the disorder strength W. 




FIG. 4: The shift of the critical energy E c vs. disorder 
strength W. (L 2 = 30 2 , <f> = 0.1, c=l). The results are plot- 
ted for the first band (Bl) and for the second band (B2). In 
the inset are represented the numerical results obtained for 
L 2 = 40 2 . 



pies. 
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